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ABSTRACT 

A hydrodynamic model for steady state, spherically-symmetric winds driven 
by young stellar clusters with an exponential stellar density distribution is pre- 
sented. Unlike in most previous calculations, the position of the singular point 
R sp , which separates the inner subsonic zone from the outer supersonic flow, is not 
associated with the star cluster edge, but calculated self-consistently. When the 
radiative losses of energy are negligible, the transition from the subsonic to the 
supersonic flow occurs always at R sp ~ 4R C , where R c is the characteristic scale 
for the stellar density distribution, irrespective of other star cluster parameters. 
This is not the case in the catastrophic cooling regime, when the temperature 
drops abruptly at a short distance from the star cluster center and the transition 
from the subsonic to the supersonic regime occurs at a much smaller distance 
from the star cluster center. The impact from the major star cluster parameters 
to the wind inner structure is thoroughly discussed. Particular attention is paid 
to the effects which radiative cooling provides to the flow. The results of the 
calculations for a set of input parameters, which lead to different hydrodynamic 
regimes, are presented and compared to the results from non-radiative ID nu- 
merical simulations and to those from calculations with a homogeneous stellar 
mass distribution. 
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Introduction 



High spatial resolution observations have modified drastically our view of star forma- 
tion in starburst and normal galaxies leading to the discovery of a large number of extremely 
massive and compact young stellar clusters with masses 10 5 M - 10 6 M , effective radii 3 pc - 
5 pc and ages less than a few times 10 6 yr - 10 7 yr. The clusters, whose masses and expected 
lifetimes coincide with those of the globul ar clusters, and which are usually named in the lit 



eratu re super star clusters (see for a review lTurnerll2009l ; Ide Grijsll2010t iPortegies Zwart et al. 



2010l ). Massive and compact star clusters are also a common feature in the nuclei of spiral 
galaxies. While the effective radii of nuclear clusters are similar to the sizes of globular and 
super star clust ers, their luminosi t ies and thus masses exceed them by up to two orders of 



magnitude (e.g. iBoker et al.l 12004 ; IWalcher et al.ll2005l ). 



Extremely high concentration (up to 10 4 — 10 5 M pc~ 2 ) of young stars implies that 
super star clusters (SSCs) are some of the most powerful feedback agents which heat up, 
shape and enrich the interstellar gas with the product of supernova explosions. The combined 
action of such clusters may lead to the formation of powerful gaseous outflows (galactic 
superwinds) which link reg ions with extreme star formation activity (starbursts) to the 



intergalactic medium (e.g . iHeckman et al.l Il990l; iMarlowe et al.l Il995t iRupke et al.l 12005 



Westmoquette et al.ll2008l see for a review I Veilleux et al.l 2005 and references therein). Thus, 
the hydrodynamic feedback from SSCs to the surrounding ISM is crucial in many respects, 
as it is the origin and the duration of the starburst event, the impact that starbursts provide 
to their host galaxies and the intergalactic medium, the evolution of the assembling galaxies 
and the feeding of supermassive black holes located within nuclear starburst regions. 



It was suggested by IChevalier fc Clegd fll985l ) (CC85) and IChevalierl fll992h that the 
kinetic energy supplied by supernovae and stellar winds is thermalized within the star cluster 
volume. This leads to a high central overpressure which drives the inserted gas away from 
the cluster in the form of a star cluster wind. If the sources of energy and mass (stars) are 
homogeneously distributed inside the star cluster volume with an outer radius Rsc and the 
radiative losses of energy are negligible, then the hydrodynamics of the outflow depend only 
on the energy and mass deposition rates, L$c and M$c, and the radius of the cluster Rsc- In 
this case the distributions of the hydrodynam i cal variables can be calcu l ated analytically (see 



Chevalier fc Clegg 



1985 



Canto et al.l 12000 ). ISilich et al.l (120031 . 12004J ). iTenorio-Tagle et al. 



f l2007h . IWiinsch et al.l J2008i boilh showed that in the case of very massive and compact 



clusters, radiative cooling becomes an important issue and developed a radiative star cluster 
wind model. 



Real clusters, however, are not homogeneous. Observations (see IPortegies Zwart et al. 
2010l and references therein) show that the stellar density drops rapidly with distance to the 
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star cluster center and that the characteristic scale for the stellar density distribution (core 
radius) is much smaller then the star cluster effective radii. This effect has been discussed by 
Rodriguez- Gonzalez et al.l (120071 ) who found an analytic non-radiative solution in the case of 
power law stellar density distrib utions and then compared it to 3D gas dynamic simulations. 
It is important to note that in IRodrfguez-Gonzalez et al.l (120071 ) solution the central gas 
density goes to infinity if the stellar density drops faster than 1/r. This implies that for steep 
stellar density distributions radi ative cooling in the central zone o f the cluster might also be a 
crucial factor. Note also that in IRodrfguez-Gonzalez et al.l (120071 ) solution the stellar density 
(and thus the energy and the mass deposition rates) abruptly drop to zero value at some 
distance from the star cluster center and it is assumed that the transition from the subsonic 



to supersonic regime occurs exactly at this point. Thus, in IRodrfguez-Gonzalez et al.l (120071 ) 
model the cut-off radius in the stellar density distribution re mains an i mport ant parameter 
although the stellar density distribution is not homogeneous. |Ji et al.l ( 120061 ) neglected the 
gravitational pull from the cluster and radiative cooling and integrated the ID stationary 
hydrodynamic equations numerically in the case of an exponential stellar density distribution. 
These results were then used to discuss the impact of non-equilibrium ionization onto the 
star cluster wind X-ray emission. 

Here we present a steady state semi-analytic hydrodynamic solution for winds driven by 
stellar clusters with an exponential stellar density distribution, which includes also radiative 
cooling. As in the above mentioned papers, we regard such a stellar distribution as more 
realistic than the formerly used homogeneous one. We thoroughly discuss boundary and 
initial conditions, which define the position of the singular point, and thus the radius where 
the flow changes the hydrodynamic regime and becomes supersonic. We also discuss the 
impact which radiative cooling provides to the flow. The paper is organized as follows: the 
input star cluster wind model is formulated in section 2. In section 3 we introduce the set of 
main equations and present them in the form suitable for numerical integration. The initial 
and boundary conditions are discussed in section 4. The method, which allows to select a 
wind solution from a branch of possible integral curves, is described in section 5. In section 
6 we first present the results of the simulations for three reference models with different 
star cluster mechanical luminosities and then discuss how other input parameters affect the 
solution. An extreme regime with catastrophic gas cooling is discussed in section 7. We 
compare our results with those obtained under the assumption that stars are homogeneously 
distributed within the cluster volume in section 8. Our major results are summarized in 
section 9. 
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The model 



We consider a young stellar cluster with constant energy and mass deposition rates Lsc 
and M$c and an exponential stellar density distribution: 

p*(r) = p* exp(-r/ J R c ), (1) 

where p*o is the central stellar density and R c is the radius of the star cluster "core". The 
total mass of the cluster is then: 



M> 



sc 



poo 

An p^r 2 exp(-r/R c )dr = 8np M Rl- 
Jo 



(2) 



Note, that the half-mass radius of the cluster in the case of exponential stellar density 
distribution is: Rhme = 2.67R C . 

It is assumed that the kinetic energy supplied by stellar winds and supernova explosions 
is completely thermalized, that the gravitational field of the cluster is negligible, and that 
the energy and mass deposition rates per unit volume, q e and q m , are in direct proportion 
to the local stellar density: 



q e (r) = q e0 exp(-r/R c ), 
q m {r) = g m0 exp(-r/i? c ), 
where the normalization constants q e0 and g m0 are: 

q e o = L sc /8nR 3 c , 
q m o = M S c/8ttR 3 c . 



(3) 
(4) 

(5) 
(6) 



3. Basic equations 

The hydrodynamic equations for the stationary, spherically symmetric flo w driven by 
stellar winds and supernova explosions if gravity is neglected are (see , for example, 



19711 : IChevalier fc Clegdll985t ICanto et al.ll2000l : ISilich et all 12004 : LJi et al||20pj ): 
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where u, P, and p are the outflow velocity, thermal pressure and density, 7 is the ratio of 
the specific heats, Q = n e riiA(T, Z) is the cooling rate and A(T, Z) is the cooling function, 
which depends on the gas temperature T and metallicity Z. Hereafter we relate the energy 
and the mass deposition rates, L sc and M SCl via the equation: 



'SC 



MscVL/2, (10) 



and assume that the adiabatic wind terminal speed, Vaoo, is constant. Thus, the model input 
parameter Vaoo defines the mass deposition rate for a given cluster mechanical luminosity 
Lsc- 

The integration of the mass conservation equation ([7]) yields: 

pur 2 = -q m0 (R c r 2 + 2R 2 c r + 2R 3 C ) exp(-r/R c ) + C. (11) 

If the density and the velocity of the flow in the star cluster center are finite, the constant 
of integration is: C = 2q m0 R^. Using this expression and taking the derivative of equation 
([H]), one can present the main equations in the form: 

du (7 - l)(q e - Q) + (7+ l)q m u 2 /2 - 2puc 2 /r 



dr p(c 2 — u 2 ) 



(12) 



dP du 

~dr = ~ pU dr ~ ^ ' ^ ^ 

2q m oR? c 



P 2 
r*u 



r 1 r 2 



1 + «: + 2 m 1 exp( - r/flc 



c 



(14) 



where c is the local speed of sound, c 2 = jP/p- Note that the central density remains finite 
and is not zero, p 7^ g cm -3 , only if the flow velocity in the star cluster center is km s _1 and 
grows linearly with radius near the center. The derivatives of the wind velocity and pressure 
in the star cluster center then are: 

^ = [(7 - l)(le -Q)~ 2fcoc 2 /3] /pc 2 , (15) 

^ = 0, (16) 
dr 

We make use of these equations in order to move away from the center and start the numerical 
integration. 



4. Initial and boundary conditions 



In the non-radiative solution the sound speed in the center is defined directly by the 
adiabatic wind terminal speed Vaoo and does not depend on the wind central density (CC85, 
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Canto et al. 2000): c A0 = [(7 - 1) /2} l / 2 V Aoo . This is not the case when radiative cooling is 
taken into consideration. Th e central gas density, p r , and the central temperature, T c , are 
then related by the equation (ISarazin &: White! 119871 ; ISilich et al.l 120041 ) : 

1/2 



Pc 



1/2 
9m0 



A(Z,T C ) ' (17) 

where c is the sound speed in the star cluster center. Equation ( 117]) shows that the 
central temperature in the radiative solution can never exceed the non-radiative value, 
Tco = P-a.Cco/'jk, where k is the Boltzmann's constant and fi a = 14/23m# is the mean 
mass per particle in the fully ionized plasma with 1 He atom per 10 H atoms. It also defines 
the thermal pressure in the star cluster center if the central temperature is known: 



P c — kn c T c — kp c T c / \i a 



(1* 



Thus, the value of the central temperature selects the unique solution from the branch of 
possible integral curves presented in Figure 1. The wind solution is the unique solution, 
which passes through the singular point, where both, the numerator and the denominator in 
equation (fl2|) . vanish and the subsonic flow in the inner zone, r < R sp , becomes supersonic 
at r > R sp . This outer boundary condition i.e. that the integral curve must pass through 
the singular point, allows one to select the value of the central temperature which leads to 
the wind solution, and also defines the singular point position. 



5. The integration procedure 

As shown in Figure 1, there are three possible types of integral curves. We will call 
them N-type, D-type and wind- type solutions, as indicated in Figure 1. The N-type solution 
is subsonic everywhere. In this case the velocity grows monotonically with distance from 
the star cluster center until it reaches a maximum value when the numerator in equation 
( I12p vanishes. For larger distances the numerator in ffl2|) becomes negative, however the 
denominator does not, and thus the velocity drops and the flow remains always subsonic. 
The D-type solution passes the sonic point at some distance from the center. At this point 
the denominator in ( Tl2|) vanishes, however the numerator does not. The integral curve then 
reaches the maximum radius and then turns back towards the center as the velocity grows. 
Note that one can easily obtain this solution, if the velocity is used instead of the radius 
as the independent variable. The wind-type solution is the only solution, which meets the 
singular point where both, the numerator and denominator in the equation of motion vanish. 

Which type of the integral curve is selected, depends on the value of the central tem- 
perature. D-type solution occurs when the central temperature exceeds that which results 
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Fig. 1. — Possible integral curves. The value of the central temperature selects the type of the 
integral curve. If the central temperature is large, the integral curve reaches the maximum 
radius r = Ro (vertical dashed line), where the denominator in equation ( Fl2l) changes sign 
and the velocity becomes supersonic. The integral curve then goes back towards the center. 
This is non physical double value solution. If the central temperature is low, the integral 
curve reaches the maximum velocity at r = (vertical dotted line), where the numerator 
in equation ( !T2|) changes the sign. The velocity then drops with distance to the star cluster 
cent er. The velocity is always subsonic in this case. This is known as a "breeze" solution 



(e.g. 



Bjorkmanlll995l ) . The wind solution is the unique solution, which passes by the singular 



point at r = R sp (vertical solid line), where both, the numerator and denominator in equation 
f fT2|) vanish and the flow becomes supersonic for larger radius. 
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into the wind-type solution. N-type solution occurs in the opposite case. This allows to 
build up a simple iteration procedure, which allows one to obtain the central temperature 
and the position of the singular point for the wind-type solution. The procedure is based 
on the bisection method and includes three integrations with three different central temper- 
atures, T max , T m i n arid Ta = (T max + T min )/2 at each iteration step. The initial values of 
T max and T min must be selected in such a way that T max is larger and T min is smaller than 
the wind-type central temperature. The central temperature, T c , for the wind-type solution 
is always between the values of T max and T min : T min < T c < T max . At the end of each 
iteration step, either T max or T m j„ is replaced with T A , what allows to narrow the interval 
for the wind central temperature and approach closer to the singular point. We usually stop 
iterations, when the position of the singular point is within an accuracy 5 < 10~ 3 , where 
5 = [(R A - Rmin) 2 + {Ra ~ Rmax) 2 } 1/2 / Ra- Rmax, R m in and R A are the radii, at which 
the denominator or numerator in the equation of motion f[T2"j) vanishes in the solutions with 
central temperatures T max , T min and T A , respectively. This procedure allows one to obtain 
the value of the wind central temperature and to localize the position of the singular point 
with high accuracy. 

However, it does not allow to pass through the singular point and obtain the runs of 
the hydrodynamic variables for r > R sp . In order to extend integral curves outside of the 
singular radius R sp and complete the solution, one must know the values and the derivatives 
of the hydrodynamical variables in the singular point. We obtain these quantities from the 
condition that at the singular point both, the numerator and the denominator in equation 
ffT2"j) vanish and thus the velocity of the flow is equal to the local speed of sound, u sp = c sp 
(see Appendix). 



6. Results and discussion 

In order to verify our model, we first compared our results with those obtained by Ji 
et al. (2006), who integrated numerically equations flZJ-flSJ) assuming that the impact from 
radiative cooling is negligible. Ji et al. (2006) obtained the singular radius R sp = 1.97 pc for 
R c = 0.48 pc, Msc = 10~ 4 M yr _1 and Vaoo = 1000 km s -1 (see section 4 and Figure 8 in 
their paper). The run of the wind velocity obtained with our code for such an input model is 
shown in Figure 1 by the solid line. It is very similar to that, obtained in the ID simulations 
(note that Ji et al. normalized their radii to the singular radius, R sp ). The radius of the 
singular point in our calculations is R sp = 1.94 pc. This implies that in the quasi-adiabati c 



case our model is in excellent (about 1.5%) agreement with the ID results of I Ji et al.l ( 120061 ). 
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Table 1: Reference models 



Model 

(i) 


Core radius 

(PC) 

(2) 


Half-mass radius 

(PC) 

(3) 


Mechanical luminosity 

(erg s" 1 ) 

(4) 


Adiabatic terminal speed 
(km s _1 ) 

(5) 


A 


1.0 


2.67 


3 x 10 40 


1000 


B 


1.0 


2.67 


3 x 10 39 


1000 


C 


1.0 


2.67 


3 x 10 41 


1000 



Figure 2 shows the results of the calculations for our three reference models presented 
in Table 1. The model clusters have the same core radius (1 pc) and adiabatic wind terminal 
speed (1000 km s^ 1 ) but different mechanical luminosities: L$c = 3 x 10 40 erg s -1 , L$c — 
3 x 10 39 erg s _1 and L$c = 3 x 10 41 erg s _1 . These mechanical luminosities correspond to 
young stellar clusters with masses about M S c = 1O 6 M , 1O 5 M and 1O 7 M , respectively. 
The distributions of velocity, temperature, density and pressure are shown in panels a, b, 
c and d, respectively. The solid, dotted and dashed lines correspond to cases A, B and C. 
In all cases the flow velocity near the center grows almost linearly with radius, passes the 
singular point at about 4 pc distance from the center and then gradually approaches the 
terminal speed value. The wind velocities are almost identical in the quasi-adiabatic models 
A and B (solid and dotted lines, respectively). 

The impact of strong radiative cooling is noticeable in the most energetic case C. In this 
case the terminal wind velocity is smaller than in the cases A and B because radiative cooling 
reduces the amount of thermal energy in the central zone, that results into a smaller wind 
energy, and thus smaller terminal speed. The impact from radiative cooling is best noticed 
in the radial profiles of temperature and thermal pressure (panels b and d, respectively). 
In the quasi-adiabatic cases A and B the temperature drops slowly with distance from the 
star cluster core showing almost un-distinguishable distributions (solid and dotted lines in 
panel b). In model C the temperature profile is different. The temperature deviates from 
the quasi-adiabatic profile significantly already at about ten times the core radius and then 
drops rapidly to the lower permitted value, of about 10 4 K, at about 35 pc from the center. 
This leads to the fast decrease in the thermal pressure (panel d) which drops more than an 
order of magnitude at this distance. Note, that thermal pressure always drops rapidly with 
distance from the star cluster center in the free wind region. Figure 2 shows also that the 
central pressure grows with the star cluster mass/power. This is also the case for the central 
density (see panel c). However, radiative cooling does not affect the density distribution 
significantly even in the most powerful case C (panel c). 

The value of Vaoo parameterizes the ratio of the star cluster mechanical luminosity 




Fig. 2. — The stationary wind solution. Panels a, b, c and d present the wind velocity, 
temperature, density and pressure, respectively. Solid, dotted and dashed lines present the 
results of the calculations for models A, B and C. Thick and thin lines in panel d display 
the thermal and ram pressure, respectively. 
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to th e mass deposition rat e. Both parameters, L$c and Msc, vary as the cluster evolves 

999) . Beside this, the flow may be mass lo aded by the gas left over 



e.g. 



Leitherer et al 



from star formation (IStevens fc Hartwelll 120031 ; ISilich et al.l |2010| ). Figure 3 presents the 



results of the calculations for three different values of this parameter: Vaoo = 1000 km s -1 , 
Vaoo = 750 km s -1 and Vaoo = 2000 km s _1 (solid, dashed and dotted lines, respectively). 
The rest of the input parameters in this set of models are the same as in the reference model 
A: L S c = 3 x 10 40 erg s -1 and R c = 1 pc. The calculated wind central temperature (see 
panel b) increases for larger Vaoo a s it is also the case in the non-radiative solution (see Canto 
et al. 2000). The calculated wind terminal speed is similar to the adiabatic value when Vaoo 
parameter is large. However, in the low velocity case, Vaoo — 750 km s _1 , the difference 
between the adiabatic and the calculated terminal speeds is noticeable. This is because the 
density in the wind, p = M(r)/4iTr 2 u(r), is larger if the selected adiabatic wind terminal 
speed is smaller (Figure 3, panel c) that leads to a stronger radiative cooling. Strong radiative 
cooling removes a fraction of the deposited energy and thus decreases the flux of mechanical 
energy and consequently the wind terminal speed. Indeed, in the quasi-adiabatic cases with 
Vaoo = 2000 km s _1 and Vaoo = 1000 km s _1 the amount of radiated energy is negligible, 
but it reaches about 14% of the deposited energy in the model with Vaoo = 750 km s _1 . 
In this case, radiative cooling leads to the rapid decrease in the temperature distribution at 
about 35 pc distance from the star cluster center as it is also the case in the reference model 
C. This results into the sharp drop of thermal pressure at the same distance (see panel d). 
Nevertheless, the position of the singular point, R sp , remains about the same, R sp ~ 4 pc. 

Finally, Figure 4 shows how the distributions of the hydrodynamical variables change 
with the star cluster core radius. This figure displays the resulting velocity, temperature, 
density and pressure profiles for clusters which have the same mechanical luminosity (L$c — 
3 x 10 40 erg s _1 ) and adiabatic wind terminal speed (Vaco = 1000 km s" 1 ), but different core 
radii: R c = 1 pc, R c = 5 pc and R c = 0.2 pc (solid, dotted and dashed lines, respectively). 
Certainly, the speed of the wind grows faster when the cluster is more compact, as it is shown 
in panel a. The central temperatures are the same in all cases. However, the temperature 
profiles are different. In the case with a larger core radius the flow is quasi-adiabatic and the 
temperature drops slowly with radius. On the other hand, the most compact model with 
R c = 0.2 pc is strongly affected by radiative cooling (dashed line in panel b). In this case 
the temperature drops rapidly and reaches the 10 4 K value at about 17 pc distance from the 
center. One can also notice the effects of strong radiative cooling on panel d, which displays 
the calculated star cluster wind thermal and ram pressure profiles. The value of the core 
radius affects also the wind central density which grows for smaller core radius, as shown in 
panel c. 




Fig. 3. — The results of the calculations for models with different Vaoo- Panels a, b, c and 
d present the wind velocity, temperature, density and pressure, respectively Solid, dotted 
and dashed lines present the results of the calculations for models with Vaoo — 1000 km s _1 , 
Vaoo = 2000 km s^ 1 and Vaoo = 750 km s _1 , respectively. Thick and thin lines in panel d 
display the thermal and ram pressure, respectively. 




Fig. 4. — The results of the calculations for models with different R c . Panels a, b, c and d 
present the wind velocity, temperature, density and pressure, respectively Solid, dotted and 
dashed lines present the results of the calculations for models with R c — 1 pc, R c = 5 pc 
and R c = 0.2 pc, respectively. Thick and thin lines in panel d display the thermal and ram 
pressure, respectively. 
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7. The catastrophic cooling regime 

The density in the wind grows with the star cluster mechanical luminosity/mass. It 
also increases if the wind is mass loaded (as the adiabatic wind terminal speed is smaller) or 
if the cluster is more compact (for a smaller R c ). In all these cases the impact of radiative 
cooling on the flow hydrodynamics becomes progressively more significant, as discussed in 
the previous section. The turn-off point in the temperature distribution moves towards the 
star cluster center and the temperature rapidly drops to the 10 4 K value. At larger radii it 
falls even to lower values because of the gas expansion. Hereafter we will assume that the 
outflow is photoionized by young massive stars and thus the wind becomes isothermal as 
soon as the temperature reaches the 10 4 K value. Outwards of this radius we replace the set 
of the main equations ffl2|) - ffT4l) with equations describing isothermal flows: 



du 2pc 2 /^r — (1 + c 2 / ^u 2 )q m u 
dr (1 — c 2 /^u 2 )pu 

2q m oB? c 

P = 2 

r*u 



fy 1 7"* 

1_(1 + i?: + 2^ )exp( ~ r/jRc 



(19) 
(20) 



where the sound speed, c 2 = jkT/fi a , is constant and the temperature is T = 10 4 K. The 
thermal pressure then is: P = pc 2 / 1 -/. However, in the numeric code we obtain thermal 
pressure by integrating the differential equation 

dP c 2 fq m p du 2p\ ^ . 



dr 7 \ u u dr r 

This allows to minimize changes in the numerical code, when the transition occurs from the 
radiative to the isothermal regime. 

In the catastrophic cooling regime the singular point detaches from its quasi-adiabatic 
position and then rapidly moves t owards the center. One can note that in this respect our 



solution is similar to that, found by lBisnovatyi-Kogan &: Blinnikovl ( Il980l ) for the accretion of 



the pre-heated gas onto a neutron star. In the later case the singular point is also located far 
away from the center, near the Bondi radius, if the accretion rate is low and the pre-heating 
of the accretion flow is small. However, when the accretion rate grows, the catastrophic 
heating regime sets in. The singular point detaches then from the Bondi radius and moves 
rapidly towards the neutron star surface. 

Figure 5 presents an example of the catastrophic cooling regime. In this case the star 
cluster core radius and the adiabatic wind terminal speed are the same as in our reference 
model A, but the mechanical luminosity of the cluster is two orders of magnitude larger: 
Lsc = 3 x 10 42 erg s -1 . The singular point moves then inside the core radius, to R sp = 0.66 pc 
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Fig. 5. — Catastrophic cooling regime. The distribution of the hydro dynamic variables in 
the model with L sc = 3 x 10 42 erg s -1 , R c = 1 pc and V Aoo = 1000 km s^ 1 . Panels a, b, c 
and d present the wind velocity, temperature, density and pressure, respectively. Solid and 
dotted lines in panel d display the thermal and the ram pressure, respectively. 
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distance from the star cluster center. The maximum flow velocity is much smaller, than in 
the quasi-adiabatic case. It reaches only about 300 km s _1 at the singular point (see panel 
a in Figure 5). The flow slows down then to about 16 km s _1 being loaded with the newly 
re-inserted matter with zero initial momentum. In this regime the last term in equation (jSj) 
dominates over the thermal pressure gradient. The velocity then increases slowly because in 
the isothermal regime the ionizing photons heat up and dump additional energy to the flow. 
The density grows in the region, where the wind slows down, and then decreases again in 
the isothermal wind region. The temperature drops abruptly to the 10 4 K value when the 
flow passes the singular point. The thermal pressure also drops just after the singular point 
and then increases in the region where the density grows up and the temperature reaches a 
constant value. 

In the case under consideration the radiative losses of energy are catastrophic indeed. 
Figure 6 shows the ratio of the energy flux through the surface with radius r: L w (r) = 
47rpu(r)r 2 (u 2 (r)/2 + H(r)), where H is the enthalpy, to the mechanical energy input rate 
inside the enclosed volume: Lscij") = LgcW — (1 + r/R c + r 2 /2R 2 C ) exp(— r/R c )]. The 
solid line in Figure 6 displays this ratio for the reference model A, whereas the dotted one 
shows the same ratio in the catastrophic cooling regime. In the catastrophic cooling regime 
L w (r)/Lsc( r ) drops very rapidly to about 10~ 3 value, whereas in the quasi-adiabatic model 
A it is close to unity at any distance from the star cluster center. This implies that in 
the extreme case, here discussed, the star cluster wind carries away only about 0.1% of the 
deposited mechanical energy. The rest is radiated away because of strong radiative cooling. 

Note, that the catastrophic cooling regime also sets in if one provides simulations for 
less energetic, but more compact clusters, and in the case of less energetic models with 
smaller adiabatic wind terminal speed parameter. In this regime the terminal wind velocity 
is small and may fall below the escape velocity. In this case a fraction of the re-inserted 
matter might remain gravitationally bound and accumulate inside the star cluster volume. 
Thus, in the cata strophic cooling re gime the gravitational pull from the cluster becomes an 



important factor (ISilich et al.ll2010l ). which should be included into the model. The impact 
that gravitational field of the cluster provides on the flow will be discussed in a future 
communication. 



8. Comparison with homogeneous model predictions 

In this section we confront the predictions from the exponential model with those, 
obtained under the assumption that stars are homogeneously distributed within the star 
cluster volume. Throughout this section we will assume that the two clusters have the same 
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Fig. 6. — The catastrophic cooling regime. The ratio of the mechanical energy flux to the 
energy input rate as a function of radius. The solid and dotted lines present the results of 
the calculations for the reference model A and the model with the same core radius and 
adiabatic wind terminal speed, but two orders of magnitude larger mechanical luminosity, 
respectively. 
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mass and their winds the same adiabatic terminal speed, but different, either exponential 
or homogeneous, stellar mass distributions. As shown in section [6] and in our prior papers, 
the distribution of the hydrodynamical variables and thus the observational manifistations 
of star cluster winds strongly depend on the characteristic space scale of the stellar mass 
distribution: the core radius, R c , in models with an exponential stellar density distribution 
and on the star cluster radius, Rsc, in models with a homogeneous stellar distribution. 
Thus, one has to link these two paramet ers in or d er to compare the models. This could 



be done in different ways. For example, I Ji et all (120061 ) compared two models assuming 



that in both cases the singular points are located at the same distance from the star cluste r 



center. One can instead use the same half- mass radius Rhm (e.g. iPortegies Zwart et al.ll2010l ). 
Specifically, here we assume equally massive clusters with different stellar mass distributions 
but with the same half-mass radii: Rhme = Rhmh, where the half-mass radius is Rhme = 
2.67R C in the exponential case and Rhmh = 0.79Rsc m models with a homogeneous mass 
distribution. The relation between the core radius R c and the star cluster radius Rsc then 
is: 

Rc « 0.3R SC . (22) 

We present three cases: our reference models A and B and an intermediate model with an 
energy input rate Lsc — 10 41 erg s _1 and assume that in each case the mass distribution may 
be either exponential, or homogeneous. The results of the calculations are presented in Figure 
7. Solid, dotted and dashed lines in Figure 7 display the results of the calculations for models 
with L S c = 3 x 10 40 erg s -1 , L S c = 10 41 erg s _1 and L S c = 3 x 10 41 erg s _1 , respectively. 
Thick and thin lines show the distributions of the hydrodynamical variables in the case 
with an exponential and a homogeneous mass distribution, respectively. One can note, that 
models with an exponential stellar mass distribution are less affected by radiative losses of 
energy. Indeed, in the calculations with homogeneous mass distribution the temperature 
and the thermal pressure already deviate significantly from the quasi-adiabatic profiles when 
the star cluster mechanical luminosity is Lsc — 10 41 erg s -1 whereas in the exponential case 
are not (compare thin and thick dotted lines in panels b and d). The mechanical energy 
input rate in the most energetic homogeneous model with Lsc = 3 x 10 41 erg s -1 exceeds the 



threshold value (see Figure 2 in iTenorio-Tagle et all 120071 ). In this case the stagnation point 



(the point where the wind velocity is km s 1 ) shifts from the center to R st = 1.9 pc and 
the shock-heated pla sma becomes thermally unstable within the central zone with r < R st 



(IWunsch et al.ll2008l ). We did not find a similar bimodal regime in calculations with an 
exponential stellar distribution. In these cases, the faster drop in density inhibits catastrophic 
cooling in the center. The two models are quite different in this respect. In the models with 
a homogeneous star distribution, the position of the singular point is fixed at r = R S c but 
the stagnation point may move from the center due to catastrophic cooling. In the models 
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Fig. 7. — The comparison of the exponential and the homogeneous model predictions. Panels 
a, b, c and d present the distributions of the wind velocity, temperature, density and thermal 
pressure, respectively. Solid, dotted and dashed lines display the results of the calculations 
for models with L sc = 3 x 10 40 erg s _1 , L sc = 10 41 erg s~ 4 and L sc = 3 x 10 41 erg s _1 . 
Thick and thin lines show the distributions of the hydrodynamical variables in the case 
with exponential and homogeneous mass distribution, respectively. It was assumed that the 
adiabatic wind terminal speed parameter is the same in all cases: Vaoo = 1000 km s _1 . 
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with an exponential stellar mass distribution it is quite the opposite: the stagnation point 
remains always at the center, whereas the singular point may detach from its quasi-adiabatic 
position and move closer towards the center. Thus the definition of the threshold mechanical 
luminosity as the mechanical luminosity above which the stagnation point moves from the 
star cluster center and the central zone becomes thermally unstable does not occur in models 
with an exponential stellar mass distribution. However, the flow may be thermally unstable 
outside of the singular point in this case. This will be thoroughly discussed in a forthcoming 
communication. 
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Fig. 8. — The distribution of the X-ray emission along the star cluster wind. The solid and 
dotted lines show the distribution of the X-ray luminosity ex (see the text) along the star 
cluster wind in the case of the exponential and the homogeneous stellar mass distribution, 
respectively. Note, that both luminosities are normalized to the total homogeneous wind 
luminosity, Lxtot = 8.1 x 10 38 erg s _1 , and that X-ray emission is slightly more concentrated 
in the case with homogeneous stellar mass distribution. 
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Diffuse X-ray emissio n has been detected from many young stellar clusters and their as- 

see the review of the 



sociated HII re gions fe.g. iMoffat et al 



resent r esults inlTownslev et al.ll201lf). 



2002 



Law fc Yusef-Zadeh 



Canto et al 



f 20001 ). 



200 



Raga et al 



(boOlUstevens fc Hartwell 



fl2003h . ISilich et all ((20051), iRockefeller et al.l ((20051) and iRodrfguez-Gonzalez et all (120071 ) 

suggested that the observed diffuse X-ray emission manifests the hot, shock-heated star 
cluster winds. The co ntribution from the hot massive stars to the observed X-ray emission 
has been discussed by lOskinoval (120051 ). The X-ray luminosity of the star cluster wind then 



is: 



Air 



Ron 



r 2 n e niA x {T, Z)dr, 



(23) 



where n e (r) and ri j(r) are the electron a nd io n number densities, Ax{Z, T) is the X-ray 
emissivity used by IStrickland fc Stevens! ((20001 ) and R out marks either the location of the 
outer wind driven shock, or the X-ray cut-off radius (the radius where the temperature in 
the wind drops below T cut ~ 5 x 10 5 K). We integrate equation (123!) numerically using the 
temperature and density profiles obtained from calculations with either exponential or ho- 
mogeneous stellar mass distribution, assuming that n e = rii = p(r)//ij, where /x, = 14/llm# 
is the ion number density. We found that the exponential model predicts a slightly smaller 
(within a factor of two) X-ray luminosity. For example, when the mechanical luminosity of 
the cluster is L$c = 3 x 10 40 erg s _1 (model A), the calculations with exponential stellar mass 
distribution predict the total 0.3 keV - 8.0 keV wind luminosity L X tot = 5.2 x 10 38 erg s _1 
whereas the homogeneous model leads to Lxtot = 8.1 x 10 38 erg s -1 . When Lsc = 10 41 erg 
s -1 , the exponential model predicts Lxtot = 6.1 x 10 39 erg s" 1 whereas the homogeneous one 
Lxtot — 1-0 x 10 40 erg s _1 . We cannot compare the X-ray luminosities in the most energetic 
case C because the central zone in the model with homogeneous stellar mass distribution 
is thermally unstable and the distributions of the hydrodynamical variables inside this zone 
cannot be obtained in the semi-analytic calculations. Figure 8 compares the distributions of 
the X-ray emission along the wind, ex = [Lxtot ~ L x (r)]/ Lxtot, in the case when the star 
cluster mechanical luminosity is Lsc = 3 x 10 40 erg s _1 . Very similar results were obtained 



in the ca 



culations where instead of Rh me = Rhmh the same singular radius, as suggested by 



Ji et al.l ((20061 ) was used for the two stellar mass distribution models. 



9. Summary 

Here we present, for the first time, a radiative semi-analytic solution for steady state, 
spherically-symmetric winds driven by stellar clusters with an exponential stellar density 
distribution. The method, here developed, improves previous calculations provided for stellar 
clusters with a given size and a homogeneous stellar density distribution and thus leads to 
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more reliable hydro dynamic predictions. It may be easily extended to clusters with other 
stellar density distributions. 

In our model, unlike in most previous calculations, the position of the singular point, 
R sp , where the transition from the subsonic to the supersonic flow occurs, is not associated 
with the star cluster edge, but calculated from the condition that the integral curve must 
pass through the singular point. When radiative losses of energy are negligible, the singular 
radius is always about R sp ~ 4R C , where R c is the star cluster core radius, irrespective of 
the other star cluster parameters. This is not the case in the catastrophic cooling regime, 
when the temperature drops abruptly at a short distance from the star cluster center and the 
transition from the subsonic to the supersonic regime occurs at the much smaller distance 
from the star cluster center. 

Radiative cooling becomes a significant factor when the cluster is very energetic /massive, 
compact, or the adiabatic wind terminal speed parameter, Vaoo = (ZLsc/Msc) 1 ^ 2 , is small. 
In the catastrophic cooling regime outflows carry away of the star cluster region only a small 
fraction of the deposited mechanical energy. The gravitationally bound, partially ionized 
nebulae may be formed then, if the photoionized gas cannot escape the gravitational well of 
the cluster. On the other hand, the low mass clusters with small energy input rates and large 
radii drive quasi-adiabatic winds. In these cases our results show an excellent agreement with 
the results of non-radiative ID numerical simulations. 

The star cluster driven wind model presented here may be applied to many problems, 
which are currently discussed in the literature. Such, as the star cluster diffuse X-ray emis- 
sion, the origin of compact HII regions, which are frequently detected around young massive 
clusters and the origin of the low-ionization line emission in the starburst driven galactic-scale 
outflows. We will address some of them in a future communication. 
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131913. GBK has been partially supported by Russian Foundation for Basic Research, grants 
08-02-00491 and 11-02-00602, the RAN Program "Origin, formation and evolution of objects 
in the Universe" and Russian Federation President Grant for Support of Leading Scientific 
Schools NSh-3458.2010.2. 



-23 - 



Appendix 

In order to obtain the flow velocity, pressure, their derivatives, and also the density and 
the temperature at the singular point, one has to use the condition that at this point the 
numerator and the denominator in equation f|T2|) vanish. The denominator in equation f|T2|) 
vanishes when the wind velocity reaches the local speed of sound and thus: u sp 
density in the singular point (see equation [H]) then is: 



c sp . The 



Psp 



1 + 



R 



R 



The second condition, that the numerator in equation (fl2~l) vanishes, then yields: 



(Al) 



ct - 2F 1 (R sp )c% + F 2 (R sp )A(T sp , Z) = 



(A2) 



where functions F± and F 2 are: 
(7-1) 
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2R r 



exp(-R sp /R c 



(A5) 



Hi = 14/ '11m h is the mean mass per ion. 



This nonlinear algebraic equation defines the temperature at the singular point, T sp , if 
R sp is known. One can present equation (IA2I) in the dimensionless form and then solve it 
numerically: 



FA = 1 - 2F 1 (R sp )c; p 2 + F 2 (R sp )A(T sp , Z)c~ p A = 0. 



(A6) 



Equation ( ]A6j) may have one, two, or have nor real roots as it is shown in Figure 9, which 
displays function FA tabulated at R sp = 4 pc radius for three different values of the star 
cluster mechanical luminosity: Lsc = 3 x 10 40 erg s _1 , Lsc = 3 x 10 41 erg s -1 and Lsc = 
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3 x 10 42 erg s 1 - solid, dotted and dashed lines, respectively. The proper solution of equation 
i6i) is selected from the condition, that segments of the integral curve obtained by the 



outward integration from the star cluster center and by the inward integration from the 
singular point match in an interior radius < Rfu < R sp . Note, that one has to obtain the 
position of the singular point, R sp , by iterations as it is described in Section 5. 



u 



sp 



Having the value of T sp , one can obtain the velocity in the singular point, which is: 
= c sp . The density in the singular point yields from equation (IA1[) . the pressure then is: 



Pap — PspC 2 sp ll- Thus, one can obtain the values of all hydrodynamic variables in the singular 
point solving the nonlinear algebraic equation (IA6j) . The value of the singular radius, R sp , 
is obtained by iterations, as explained in section [5j 
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Fig. 9. — Different roots of equation flA6j) . Solid, dashed and dotted lines display function 
FA tabulated at R sp = 4 pc for three different mechanical luminosities: L$c = 3 x 10 40 erg 



>5C 



3 x 10 41 erg s 1 and L$c = 3 x 10 42 erg s 1 , respectively. 



In order to obtain the derivative of velocity in the singular point, one can use the 
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L'Hopital's rule. The derivatives of numerator and denominator of equation ffl2|) over radius 
are: 

dN _ dN dN du dN dc 2 dN dT dN dp _ du dN 
dr dr du dr dc 2 dr dT dr dp dr 5 dr 6 dr ' 



dD dD du dD dc 2 



( +1) + ^ m ^\ + ^ 
\dr p J r 



dr du dr dc 2 dr 
where functions N, D, dN/dr, F§ and Fq are: 

N(r,u,p,c,T) = ( 1 -l)( qe -Q)-Aq m0 c 2 {R c /rf + 

g m [( 7 + l)u 2 /2 + Ac\R c /rf{l + r/R c + (r/R c ) 2 /2)} 
D(u, c) = c 2 — u 2 
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(A14) 



One can obtain then the derivative of the wind velocity (and thus the derivative of the 
thermal pressure) at the singular point substituting relations (1A7I) and (1A8j) into equation 
(fl2"|) and keeping in mind that at the singular point u sp = c sp . This leads to a quadratic 
algebraic equation: 



du\ ^ du „ 
— -2F 7 — + F 8 = 
dr J dr 



(A15) 
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where functions Fj and F$ are: 

F 7 = [2pc 2 /r - (7 + l)cq m - F 5 ] /2( 7 + l)cp . 
F 8 = [F 6 + 0jV/0r]/[( 7 + l)cp]. 



(A16) 
(A17) 



The root of equation ( IA15j) . which leads to the positive derivative du/dr at the singular 
point, is used in the calculations. 
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